Exact Dirac equation calculation of ionization induced by 

ultrarelativistic heavy ions 

Q\ ■ A. J. Baltz 

Physics Department, Brookhaven National Laboratory, Upton, New York 11973 

c. 

^ ■ (February 2, 2008) 

00 
(N 

^ ■ Abstract 

The time-dependent Dirac equation can be solved exactly for ionization in- 
^ I duced by ultrarelativistic heavy ion collisions. Ionization calculations are 

'^ ' carried out in such a framework for a number of representative ion-ion pairs. 

i-G I For each ion-ion pair, the computed cross section consists of two terms, a con- 

stant energy independent term and a term whose coefficient is In 7. Scaled 
^ \ values of both terms are found to decrease with increasing Z of the nucleus 

(^ ' that is ionized. 

O ■ PACS: 34.90. -Kq, 25.75.-q 

G^ ■ 
o 

• 1— I 

>>: 

Ph. 



I. INTRODUCTION 



In a recent work [|[| an exact semiclassical Dirac equation calculation of ionization prob- 
'k> \ abilities was presented. The case considered was colliding Pb + Pb at ultrarelativistic 

c^ \ energies. A single electron was taken to be bound to one nucleus with the other nucleus 

completely stripped. The probability that the electron would be ionized in the collision was 
calculated as a function of impact parameter, but no cross sections were presented. In this 
paper the approach of Ref. [|I| is extended in order to calculate ionization cross sections 
for a number of representative cases of collisions involving Pb, Zr, Ca, Ne and H ions. In 
Section II the exact semiclassical method is reviewed and calculations of impact parameter 



dependent ionization probabilities are presented. In Section III the results of the probability 
calculations are used to construct cross sections for various ion-ion collision combinations in 
the form 

a = A\n-f + B (1) 



where A and B are constants for a given ion-ion pair and 7 (= l/\^l--v^) is the relativistic 
factor one of the ions seen from the rest frame of the other. Comparisons are made with 
previous ionization calculations. In Section IV the example of Pb + Pb at the AGS, CERN 
SPS, and RHIC is worked out and the CERN SPS case is compared with data. 

II. IMPACT PARAMETER DEPENDENT PROBABILITIES 

If one works in the appropriate gauge [Q, then the Coulomb potential produced by an 
ultrarelativistic particle (such as a heavy ion) in uniform motion can be expressed in the 
following form Q] 

V{p, z, t) = -6{z - t)aZp{l - a,) In ^^ ~/^ . (2) 

b is the impact parameter, perpendicular to the 2;-axis along which the ion travels, p, z, and 
t are the coordinates of the potential relative to a fixed target (or ion), a^ is the Dirac matrix, 
a is the fine structure constant, with Zp and v the charge and velocity of the moving ion. 
This is the physically relevant ultrarelativistic potential since it was obtained by ignoring 
terms in (b — p)/7^ 0] 0. As will be shown in Section II, when b becomes large enough 
that the expression Eq.(2) is inaccurate, we match onto a Weizsacker- Williams expression 
which is valid for large h. Note that the 6^ in the denominator of the logarithm in Eq.(2) is 
removable by a gauge transformation, and we retain the option of keeping or removing it as 
convenient. 

It was shown in Ref. [Q that the 5 function allows the Dirac equation to be solved exactly 
at the point of interaction, z = t. Exact amplitudes then take the form 



aUt = cx)) = %+ / dte'^^f-^^^'{(j)f\6{z - t)(l - a,) 



iaZp In (b-p)2 



x(e-'"^^^°^^-^^ -1)10,) (3) 

where j is the initial state and / the final state. This ampltude is in the same form as the 
perturbation theory amplitude, but with an effective potential to represent all the higher 
order effects exactly, 

V{p,z,t) = -t6{z-t){l - a,)(e-^"^-"^(^-^)' - 1), (4) 

in place of the potential of Eq.(2). 

Since an exact solution must be unitary, the ionization probability (the sum of probabil- 
ities of excitation from the single bound electron to particular continuum states) is equal to 
the deficit of the final bound state electron population 

T.Pib) = l-Y.Pib) (5) 

ion bound 

The sum of bound state probabilities includes the probability that the electron remains in 
the ground state plus the sum of probabilities that it ends up in an excited bound state. 
From Eq. (3) one may obtain in simple form the exact survival probability of an initial state 

PAb) = 1(0.1(1 - a.)e-^-"^(^-P)^|0,)p. (6) 

In principle the ionization probability can be calculated without reference to final continuum 
states. In practice ionization will be calculated both as a sum of continuum probabilities as 
well as the deficit of bound state probabilities. 

Table I shows the results of ionization calculations for an ultrarelativistic Pb + Pb 
collision. One of the Pb ions has an electron initially in the ground state. The other is 
completely ionized. Final state probabilities for the electron are calculated as a function of 
impact parameter b. Calculations have been carried out with a logarithmic spacing in values 
of b, with the range of b chosen to go from constant probabiity of ionization at the low end 
to 1/6^ behavior at the high end. The last column which is the sum of final bound state 



(column 3) and final continuum state (column 4) population exhibits a small deficit from 
unity, presumably mostly from the truncation of the energy sum over excited continuum 
states or secondarily from the truncation of the bound state sum in the calculations. 
Tables II and III show corresponding calculations for Ca + Ca, and Ne + Ne. 

III. CROSS SECTIONS 

The actual cross section comes from the impact parameter integral 

ai^ = 271 J Pib) b db. (7) 

As was exemplified in Tables I-III, for each ion-ion case calculation of probabilities was 
performed at ten impact parameters, in a scheme of sequential doubling. The points ran 
from an impact parameter small enough that the probabilities were approximately constant 
with b, to an impact parameter large enough that the probabilities had started to fall off as 
1/6^. The part of the integral, Eq.(7), over this region from the first to the tenth impact 
parameter takes the form of sum of nine integrals on a logarithmic scale 



ai_io = E 2vr / '^' P{b) b db 

i=l,9 •'''- 



^ E 2- < mh' > T" y • (8) 

Approximating < P{b)b^ > over each interval by the average of the end points we have (for 
6i+i = 2bi) 

ai_io = 7rln2 E iPih)b^ + P(&m)&m)- (9) 

j=l,9 

Since the probability goes to a constant at the lowest impact parameter, bi, the contri- 
bution to Eq. (7) from zero up to 6i is given simply by 

ao_i = nP{b,)bl (10) 

We now need the contribution from the last point computed, 6io, out to where P{b) cuts 
off. We match a Weizsacker- Williams type calculation to the exact calculation at this 6io 



in order to determine the normalization for the calculation of probabilities at larger impact 
parameters and to make the high end cutoff in b. How the calculation in the delta function 
gauge and the calculation in the Weizsacker- Williams formulation are equivalent at this 
matching impact parameter 610 is presented in Appendix A. 

The Weizsacker- Williams expression for a transition probability at a given impact pa- 
rameter is of the form 

Pwwib) = r P{u,)^Kl{-)du:, (11) 

Jeb 7 7 

with Eb the ground state electron binding energy (a positive number here). 
If 6cj << 7 then 

^?(^) - i. (12) 



7 uj'^lP ' 



and 



1 r°° 
Pwwib) = -J^ P{uj)du. (13) 

At the matching impact parameter, Eqns. (12) and (13) are accurate up to the point where 
the energy starts to cut off. Thus /^ P{uj)duj may be simply determined 

/•oo 

/ P{iu)diu = bl,Piho). (14) 

JEb 

Next recall that to high degree of accuracy 

4rA-?(^).di,=/„(^)^ (15) 

7 •'bm 7 ^OlO 

Then from Eqns. (7), (11), (14), and (15) the contribution to the cross section for impact 
parameters greater than 610 is 



/ .68I7 f°° f°° ^ 

o"io-oo = 27r ln(— ) / P{uj)duj — / P{uj)\nujduj 

\ Oin JEb JEb I 



= 27r6?oP(feio)(ln(^)- < Inc. >j. (16) 

< lnt<; > can be evaluated from the empirical observation that at 610, -P(^) goes as l/u"' 
with n ^ 3.8. One obtains 



<lna; >= 111^5 + (17) 

(n-1) 

One now has the full ionization cross section 

(T = (To-l + O-l-lO + 0"lO-oo (18) 

or in the usual form 

a = A \n-f + B. (19) 

with 

A = 27rbloP{bio). (20) 



and 



B = nP{bi)bl + 7rln2 5] {P{b,)b^ + P{b,+i)b^ 



i+i) 



i=l,9 



+2vrP(6io)&?o (in(^) -\nEs- j^^ . (21) 

The A In 7 term is entirely from the non-perturbative, large impact parameter region 
and gives the beam energy dependence arising from the impact parameter cutoff at 6 ^ 7/u; 
(see Eqns. (15) and (16)). Despite the form of Eq.(20) A does not really depend on the 
matching impact parameter 610 since 610 is in the region where P{b) ~ 1/6^. The B term 
is independent of beam energy and contains non-perturbative components from the smaller 
impact parameters. 

Table IV shows the results of calculations of the ionization cross section components A 
and B for symmetric ion-ion pairs. There is good agreement between the cross sections 
calculated by subtracting the bound state probabilities from unity (first rows) or calculated 
by summing continuum electron final states (second rows). The agreement with the Anholt 
and Becker calculations ^ in the literature is good for the lighter species for both A and 
B. However with increasing mass of the ions the perturbative energy dependent term A 
decreases in the present calculations whereas it increases in the Anholt and Becker calcula- 
tions. The greatest discrepancy is for Pb + Pb, with Anholt and Becker being about 60% 



higher. Perhaps this discrepancy is due to the fact that Anholt and Becker use approximate 
relativistic bound state wave functions and the present calculations utilize exact Dirac wave 
functions for the bound states. Surprisingly, it is the term B (which has the non-perturbative 
component) where agreement is relatively good between Anholt and Becker and the present 
calculations of Table IV. 

Table V shows results of the calculation of B (multiplied by Z|/Z^) for a number of 
representative non-symmetric ion- ion pairs. (Since A is perturbative, scaling as Zf, its 
value can be taken from Table IV for the various pairs here.) Note that if one goes to the 
perturbative limit for a Pb target to be ionized (H + Pb or Ne +Pb) then the scaled B values 
(17,090, 17,030) are some 30% higher than the necessarily perturbative Anholt and Becker 
value of 13,000. The good agreement of Anholt and Becker with the present calculations for 
Pb + Pb B seen in Table IV is thus somewhat fortuitous. 

IV. AN EXAMPLE 

Table VI presents calculated ionization cross sections for Pb + Pb at the AGS, the CERN 
SPS, and RHIC CERN SPS data of Pb with a single electron impinging on a Au target has 
recently been published by Krause et al. 0. Their measured cross section of 42,000 barns 
is significantly smaller than the Anholt and Becker calculation (which includes screening in 
the target Au) of about 63,600 barns. The result of the present Pb + Pb calculation (55,800 
to 58,200 barns) does not include screening and should be compared with the corresponding 
no-screening calculation of Anholt and Becker (83,700 barns). What was essentially the 
present Pb + Pb result was privately communicated to Krause et al., and they seem to have 
then assumed that if screening were to be included, then the present calculation should be 
scaled by the ratio of the Anholt and Becker screened to unscreened results. They comment 
in their paper, "With screening included [Q and scaled to a Au target, the Baltz value agrees 
with the ai measured in the ionization experiment (4.2 x 10^ b)." 
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APPENDIX A: RELATIONSHIP OF THE S FUNCTION POTENTIAL TO THE 
WEIZSACKER- WILLIAMS FORMULATION 

The Weizsacker- Williams cross section for a process induced by a heavy ion projectile 
of charge Zp,(T^^(cj) is expressed in terms of the photoelectric cross section aph{uj) for the 
same process. 



oo 



<yZ„ u , , ^^o,huj , 



<y^^{u) = 27r / :-^-ap^{u)Kl{—) b db (Al) 

Now since the photoelectric cross section for a process is given by 

a,,{u;) = 1 / d^rfC^ ■ eiJoe'n', (A2) 

the Weizsacker- Williams amplitude (apart from an arbitrary constant phase) for the process 
is 

2aZr, ^^ ,bu , 



7 7 

Consider the 6 function gauge 



^Ki{—) I d^r^}a ■ bV-oe*"' (A3) 



V{p, z, t) = -5{z - t)aZ{l - a,) In ^^-^- (A4) 



Its multipole expansion is 



V{p, z, t) = aZ{l - a,)6{z - t) 

- In ^ p>b 

.^-^ 2cosm(/) 



m>0 "^ 



X 



P\"' U 

b) P^' 

p>b (A5) 



For b» p 



■Jj\m 

+ ' 



V(p, z, t) = 6(z - t)aZ(l - a,)2y cos 0. (A6) 





One may make a gauge transformation on the wave function 

^ = e-^'^('"'*V (A7) 

where 

Xir,t) = -2e{t-z)aZj-cos(j). (A8) 

This leads to added gauge terms in the transformed potential 

- M^ _ a ■ Vx(r, t) = 26(z - t)(l - aJaZ^ cos</. + 26(1 - z)aZ^^. (A9) 

at b b 

(This is the same transformation as previously carried out without the restriction b » p 
to go to the light cone gauge 0.) Here we obtain the light cone gauge potential for 6 >> p 

\/(p,z,t) = 20(t-;z)aZ^— , (AlO) 

and we then obtain the perturbative amplitude in the light cone gauge 

acone = ^j^ dtj d^rjf^-Woe''^'. (All) 

Integrate over t 

f £r'ip}a-h'ilJoe''^'. (A12) 



'^cone 



cub 
Now consider a^w For ■j » bu 



and Eq. (A3) becomes 



M-) = f . (A13) 

7 buj 



a^^ = ^^ / d^rij*fa ■ Voe^^^ (A14) 

tub J ■' 



Thus if one transforms from the delta function gauge to the light cone gauge the am- 
plitude in that light cone gauge is found to be equal to the Weizsacker- Williams amplitude 
(within an arbitrary constant phase) as long as 6 >> p and 7 >> boo. 

10 
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TABLES 



TABLE L Ionization and Unitarity: Probabilities for Pb + Pb 



b(fm) 


e 


EtodC 


2-icont ^ 


Ee- 


62.5 


.4344 


.5337 


.4399 


.9736 


125 


.4467 


.5474 


.4228 


.9703 


250 


.4884 


.5920 


.3788 


.9708 


500 


.5820 


.6828 


.2907 


.9735 


1000 


.7303 


.8165 


.1691 


.9856 


2000 


.8899 


.9447 


.0526 


.9973 


4000 


.97056 


.98986 


.00987 


.99972 


8000 


.99270 


.99777 


.00217 


.99994 


16,000 


.998178 


.999547 


.000529 


.999986 


32,000 


.999545 


.999865 


.000131 


.999996 



TABLE IL Ionization and Unitarity: Probabilities for Ca + Ca 



b(fm) 


^gr 


Ebwe 


l^cont ^ 


Ee- 


250 


.95335 


.96291 


.03548 


.99839 


500 


.95351 


.96297 


.03516 


.99812 


1000 


.95657 


.96562 


.03278 


.99841 


2000 


.96355 


.97170 


.02707 


.99878 


4000 


.97583 


.98293 


.01635 


.99927 


8000 


.99028 


.99520 


.00462 


.99982 


16,000 


.99760 


.99924 


.00074 


.99998 


32,000 


.999419 


.999830 


.000165 


.999996 


64,000 


.9998559 


.9999585 


.0000404 


.9999989 


128,000 


.9999640 


.9999897 


.0000100 


.9999997 
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TABLE III. Ionization and Unitarity: Probabilities for Ne + Ne 



b(fm) 


Cgj. 


Efen^e 


Z^cont ^ 


Ee- 


500 


.98814 


.99058 


.00901 


.99959 


1000 


.98816 


.99058 


.00895 


.99952 


2000 


.98894 


.99125 


.00830 


.99960 


4000 


.99070 


.99278 


.00692 


.99970 


8000 


.99383 


.99564 


.00418 


.99982 


16,000 


.99753 


.99878 


.00117 


.99996 


32,000 


.999392 


.999809 


.000186 


.999995 


64,000 


.9998534 


.9999572 


.0000416 


.9999989 


128,000 


.9999636 


.9999895 


.0000102 


.9999997 


256,000 


.99999092 


.99999740 


.00000254 


.99999993 



TABLE IV. Calculated Ionization Cross Sections Expressed 


in the Form Alnj + B (in 


barns) 




Pb + Pb 


Zr + Zr 


Ca + Ca 


Ne + Ne 


H + H 


1 - T,bnd e" 


8680 


10,240 


10,620 


10,730 


10,770 


^ l^cont 6 


8450 


9970 


10,340 


10,440 


10,480 


Anholt & Becker g 


13,800 


11,600 


10,800 


10,600 


10,540 


1 — Y^bnd s 


14,190 


28,450 


38,010 


46,080 


71,090 


B Econt e" 


12,920 


27,110 


36,530 


44,430 


68,780 


Anholt & Becker 


13,000 


27,800 


37,400 


45,400 


70,000 
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TABLE V. Calculated values of the scaled quantity (Zl/Z^)^ for non-symmetric combinations 
of colliding particles. The second nucleus (Z2) is taken to be the one with the single electron to 
be ionized. Since Anholt and Becker cross sections without screening are completely perturba- 
tive, their values of of B also can be taken from Table IV, and are repeated here for convenient 
comparison. 





H + Ne 


H + Ca 


Ca + H 


H + Zr 


H + Pb 


Pb + H 


1 - T,bnd e 


46,150 


38,270 


70,820 


29,440 


17,090 


67,550 


2^cont ^ 


44,490 


36,790 


68,520 


28,070 


15,680 


65,330 


Anholt & Becker @ 


45,400 


37,400 


70,000 


27,800 


13,000 


70,000 




Pb +Ne 


Ne + Pb 


Pb + Ca 


Ca + Pb 


Pb + Zr 


Zr + Pb 


1 - Efend e- 


42,560 


17,030 


34,720 


16,870 


26,010 


16,250 


2^cont ^ 


41,000 


15,690 


33,330 


15,530 


24,730 


14,930 


Anholt & Becker 


45,400 


13,000 


37,400 


13,000 


27,800 


13,000 



TABLE VI. Example: Calculated Ionization Cross Sections For Pb + Pb (in barns) 



1-E 



bnd' 



E 



cont '-' 

Anholt & Becker @] 

Anholt & Becker (with screening) 

Krause et al. Pb + Au data iHl 



S 7 = 11.3 


CERN 7 = 


160 


RHIC 7 = 23, 000 


35,200 


58,200 




101,400 


33,400 


55,800 




97,800 


46,700 


83,700 
63,600 
42,000 




151,600 
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